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Dans ce chapitre, on désigne par K le corps des réels R ou celui des complexes €, et 
par E un espace vectoriel normé de dimension finie sur X. On note aussi (L(E),+,.,0) 
l’algèbre de Banach des endomorphismes de Æ. Nous allons nous intéresser aux équations 
différentielles linéaires à coefficients constants dont le second membre est une fonction 
exponentielle polynôme. Ce chapitre constitue une application directe du Chapitre sur la 
réduction des endomorphismes. 


1 Exponentielle d’un opérateur 


. ce mn . . D. ‘a 
Le rayon de la série entière Ÿ[ # est infini de sorte que cette série converge sur tout R 
n>0 
: re lu” ’ js u” 
(vers e*). Si u € L(E), la série ÿ! F- convergera, et par conséquent la série ÿ° # sera 
n2>0 n>0 
normalement convergente sur toute boule fermée B de centre O de L(E). On définit ainsi 


un opérateur noté e“ et appelé l’exponentielle de x : 


Théorème 1 On a 

(1) e° = Id. 

(2) Vu,veL(E) -uv=vusentr— ete, 
En particulier e“ est inversible et (e“) = e 
(3) Vu,ve L(E) uv = vu ue = etu.(4) Vue L(E) {|le“| < elul. 


—U 


Preuve : Seul (2) présente une difficulté. On a 


N N 
nu . u" uP : uv? 
ee” = lim ) — ) — | — lim An avec ÀAn — ) ——, 
N—+00 n! p! nlp! 
n—= 


N—+o 
O<n<N 
OPEN 
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N k N n 
u + ; uv? 
eurv — ) ur = RL Bn avec By = ) ) ré 
k=0 ! —+00 DER np: 


Si l’on représente dans le plan les points de coordonnées entières (n,p) sur lesquels on 
somme, on constate que An correspond à un carré et que Bn correspond à un triangle 
contenu dans ce carré. Par suite 


lol LS Cul + Het) 
An — Bnl < > > nl > k! 


k=N+1n+p=k k=N+1 


se ; DS n 
et cette dernière somme tend vers 0 quand N tend vers +oo, puisque la série D 


n>0 
converge pour tout x € R.m 
Théorème 2 Soit ue L(E). L'application 
J: (R,+) — (L(E),o) 
t + ett 
est un homomorphisme de groupes dérivable sur R, et tel que f'(t) = ueït = ettu. 
Preuve : 
ft+h)-f() HR QUE CRUE T" “er (hu)" 
a — ——— — € —————>— | = — 
h k k h n! 
n>1 
hu hr2u 
= tu L= tu tu 
= EN ni = tu + th NT ZE À 
n>1 2 
—2u" 


est bornée quand |h| < 1 (on à en effet la majoration bu” 


conséquent lim ICO elly. M 


Théorème 3 Soient M : R — L(E) etx : R — E deux fonctions dérivables (resp. de 
classe C!). Alors l'application 


g: R — E 
t + Mi(t)x(t) 


est dérivable (resp. de classe C1) et g'(t) = M'(t)x(t) + M(t)x'(t). 


Preuve : La démonstration est celle de la dérivabilité d’un produit. Pour t fixé, il faut 
montrer que 


g(t+h)—39() 


A = : — (M'(t)z(t) + M (+) x (6) 


tend vers 0 quand h tend vers 0. On a 


A 


RIM + a+ h) M 9x (MO x 1 G 0) 
= MG+R) D En 
_ M (+1) [ED 0 UFR 20 | + [RE NO tin M (t) D 


Majorons |A] par la somme des trois normes, et donnons-nous un réel £ € R*. La 
fonction k + | M (t+h)| est bornée au voisinage de 0 car M est continue en t. Il existe 
donc a € R* et 7, € R* tels que |[M(t+h)| < a pour tout À vérifiant |A] < 71. Par 
hypothèses, on à ER) #6) — x (t) | < € dès que |A] < M, en va — M! C] < € 
dès que |hl] < 73, et [M (t+h) — M(t)| < € dès que [[h]| < 74. Finalement, on aura bien 
l'implication 


VEERT = Inf(m,M0304) hl <n= IA < (a+: (+12 (O1) em 


2 Systèmes linéaires d’ordre 1 sans second membre 


2.1 Théorème fondamental 


Un système différentiel linéaire d’ordre 1 sans second membre est la donnée d’une 
équation 


'(t)=ux(t) (4) 
où u € L(E) est fixé. L’inconnue x (t) est une fonction dérivable t + x (t) de R dans 
E. Toute fonction dérivable t + x (t) de R dans Æ vérifiant (A) est appelée solution du 


système différentiel (H). Si l’on rapporte E à une base e = (e1,.….,en) et si l’on note 
(&i;) la matrice de u dans e, le système différentiel (A) s'écrit 


æ' (€) = Mzx(t) 


où æ(t) = t(xi(t),.….,29 (t)) désigne le vecteur colonne des coordonnées de x (+) dans la 
base e. In extenso, le système (A) s'écrit 


Le problème de Cauchy relatif à l'équation différentielle (H) et aux conditions initiales 
(to, xo) € K x E est celui de la recherche des solutions du système (A) vérifiant x (to) = o. 


Théorème 4 1) Le système (H) : x!(t) = ux(t) admet la fonction x (t) = et-to)u (x) 
pour unique solution satisfaisant la condition initiale x (to) = to. 
2) L'ensemble S des solutions de (H) est un K-espace vectoriel de dimension n. 


Preuve : 1) La fonction æ(t) = elt-t0)(x5) est clairement solution de (H) et vérifie 
x(to) = xo. Si x(t) est une autre solution de (H), posons g(t) = e t““x(t). On a 
g'(t) = —ue tüx(t) +e tx! (t) = 0, de sorte que g(t) soit une fonction constante. En 
posant g(t)=ce E, on obtient x (t) = el“ (c). La condition initiale x (to) = xo implique 
c= et (xç), de sorte que x (t) = et" (e—tot (x)) = ebé-to)u (x5). 

2) L'ensemble S est clairement un K-espace vectoriel et l’application 


wo: E — $S 
to + æ(t) = e“*(xo) 


est un isomorphisme d’espaces vectoriels. En effet, est linéaire, surjective puisque toute 
solution x dans S s'écrit x (t) = e!" (x) avec x0 = x (0) d’après le 1), et injective car 


VER zx(t)=et"(xo) =0 = x9 = 0. 


2.2 Structure des solutions complexes 


Soit u € L(E). Notons À, …, À, les valeurs propres de u, puis @1, …, a, les multipicités 
correspondantes de ces valeurs propres dans le polynôme caractéristique x,, de u, et enfin 
N (x) = ker (u — À;1d)" le sous-espace caractéristique associé à À;. On sait que 


CL) dim N (x) — di, 
(2) E=N(xM)8.….8N(h), 
(3) AN) est stable par u. 


Tout vecteur x de E se décompose en x = y1 +. + y avec y; € N (À), et donc 


P P 
ett (x) — Je. (y) = de rate (y) 

i=1 i=1 
p œil k 

_ DCS #k Cu (y) car etXild = tXIq 
i=1 k=0 Â 
P 

— DD A0 et 
i=1 


œil EN 

où PR) S GX Ia) (y:)t* est un polynôme en t à coefficients dans E. On peut 
k=—0 

énoncer 


Théorème 5 - Structure des solutions complexes - 
p 

Les solutions de (H) : x'(t) = ux(t) sont de la forme x(t) = D P;(t)et\ où les À; 
i=1 


désignent les valeurs propres de u et où les P;(t) sont des polynômes de E[t] de degrés 
deg P; < @;, @ représentant la multipicité de À; dans le polynôme caractéristique de u. 


Théorème 6 Si toutes les valeurs propres de u sont simples, les solutions de (H) sont 


p 
exactement les fonctions x (t) = D yiet\ où y; décrit N (X). 
i=1 


2.3 Calcul pratique des solutions complexes 


p 
On remplace x (t) par l'expression Ÿ P, (t) et dans l'équation (H), pour obtenir 


i=1 
P P 
D(PFORARO = TS UP te" 
i=1 em es —, 
EN(X) EN(X) 


Les appartenance à N (;) se vérifient en remarquant que tous les coefficients des P; 
sont dans N(;), N'(À;) étant stable par u et donc par u — À;1d. La somme directe 
E=N(X)®..@Nt(x) permet d'écrire 


P; (4) + AP, (6) =u(P, (t)) 


pour tout indice 4, et cela déterminera parfaitement les polynômes P;. Cette méthode de 
résolution est dite ” par identification”. 


2.4 Solutions réelles 


Si (A) est un système différentiel réel, la matrice de u est à coefficient réels. Le polynôme 
caractéristique de « l’est aussi et ses racines sont soit réelles, soit conjuguées. Notons 
A1,.…., À9 les valeurs propres réelles de u, a1,.…., a, leurs multipicités, p,.…., 4, 1,7 les 
valeurs propres complexes non réelles et 71,...,7, les multipicités de w1,.….,u,. l'application 
u : R° — KR’ est aussi une application de C” dans C”, et l’on connaît parfaitement les 
solutions complexes de (H). Pour obtenir toutes les solution réelles de (H), le plus simple 
est d'utiliser le 


Lemme 1 La fonction x (t) est une solution réelle de (H) : x’ (t) = ux (t) si, et seulement 
si, c’est la partie réelle d’une solution complexe. 


Preuve : Toute solution réelle est bien la partie réelle d’elle-même. Réciproquement, si 
æ(t) = y(t) +iz(t) est solution complexe de (H), alors y (t) + iz/ (t) = uy(t) + iuz (t) 
entraîne y (t) = uy(t), et y(t) sera bien solution de (H). m 

Une autre méthode possible consiste à donner à priori la forme des solutions, puis à ré- 
injecter l'expression dans l’équation différentielle. Cette méthode, dite ’méthode par 
identification”, est basée sur le Théorème suivant : 


Théorème 7 - Structure des solutions réelles - 
Les solutions réelles de l'équation (H) : x’ (t) = ux(t) sont de la forme 


P 


q q 
2 (6) = TP ef + DO (He + SG 
i=1 i=1 


i=1 


où P;(t) € R'[t] sont des polynômes de degrés deg P; < œ, et où Qu(t) € C”[t] avec 
deg Qi < Yi. 


Preuve : Tout vecteur x € E se décompose de façon unique en 
= sdb hear due (1) 
avec % E NX), € Nu) et z € N(%;) donnés par le Théorème 5. Si x € R”, 
TS 2 7 en Gt nn (2) 


Les équations (1), (2) et le Lemme ci-dessous entraînent y; = 7, et z; = 7! pour tout 4, 
soit y; € R et 


T=Yte EYp ++ +2 +7 +. +. 


Toute solution réelle s’écrira donc 


P 


q 
ett (x) — ÿ et" (y) + D» et (2) + D? eu (z;). 
i=1 


i=1 


En refaisant le même calcul qu’au Théorème 5, on trouve 


P q q 
st) = NP (ter + NT Q; (te + SUR; (t) ef: 
i=1 i=1 i=1 


avec 
Vi LL CLDE V1 LC TPE 
QD = DUR €) où RD = D CRD €) 2 GT 
k=0 k=0 


de la forme désirée. = 


Lemme 2 On a N{u) = N(f%) pour toute valeur propre a de u. 


Preuve : Les racines y et x ont même multiplicité a dans x,,, donc 


xEeN(u)&(u-uld) (x) =0£&(u-x1d) (F)=087TeEeN(x) 


ce qui entraîne MN (u) C N(x). Réciproquement, tout y € N'(%) s'écrit y = T et les 
équivalence ci-dessus entraînent x € N (u), soit y € N(u).m 


2.5 Etude des systèmes différentiels 2 x 2 


Etudions les trajectoires des solutions de l’équation 


(H) zx'(t) = Mxt(t) 


où la matrice M — ( : : ) est à coefficients réels. 


a) Si M possède deux valeurs propres distinctes À et u (À < y), alors M est diagonalisable 
et les solutions de (H) sont les fonctions x (t) = ceux + coeltu3 où (ui, u3) désigne une 
base de vecteurs propres de M. La trajectoire des solutions est l’arc paramétré 


{ 21 (4) = ae’ 


x (t) = coelt 


et l’on obtient l’une des allures de courbes suivantes : 


€ = 0, &@ >0et nu <0 a =0,@>0etu>0 C0 £0et À<yu<0Û 


Si c1c2 £ 0, l’origine est stable si À < y: < 0 (cela signifie que l’on est amené vers l’origine 
lorque le paramètre t croît), et instable si 0 < À < u. Dans les deux cas, les trajectoires 
sont les mêmes mais sont parcourues en sens inverse. 

b) Si M possède deux valeurs propres réelles confondues, on envisage deux cas : 


e Si M est diagonalisable, les solutions de (H) sont les fonctions x (4) = ce tui + ce tu 
où (4j, u2) est une base de vecteurs propres de M. La trajectoire est une demi-droite de 
vecteur directeur cui + cou3. 

e Si M n’est pas diagonalisable, alors dim Ker (M — Ad) = 1 et 


E = Ker (M — Xd)° 3 Ker (M — Xd) = E(X). 


Si (ui, u2) est une base de E telle que ui € E (À) et u3 € E\E (À), une solution de (H) 
s'écrit 


a (t) = et (0) = Met Ad) (ont + cou) = ect + e*co (m8 + t(u — XId) ü) 


Comme (u — X1d)u3 € E(À) entraîne l’existence de a € R tel que (u — A1d)u3 = au, 
on obtient x (4) = (c1 + cat) etui +e"œus et les trajectoires des solutions admettent la 
paramétrisation 


Les allures des trajectoires seront : 


À>0 À<0 


c) Si M possède deux valeurs propres complexes non réelles À = à + if et À = a — ip, 


la forme générale des solutions réelles est x (4) = ye" + %et où y = a + 1b, a,b € R?° et 


y € E(À). Comme7e E (À) et C?=E(A@DE (X), le couple (y,7) est une base de C? et 


(a,5) = (Rey, Im y) = (52 13) 


sera une base de C? formée de deux vecteurs réels. Dans cette base (a, b) : 
a (t) = (a +ib) e%eiPt + (a — ib)e%e 8 = 2% (a cos Bt — bsin Ft) 
et la représentation paramétrique d’une trajectoire t + x (t) sera de la forme 


1 (t) = 2e% cos Bt 
xo (t) = —2e% sin Bt. 


Le diagramme est stable si à < 0, instable si & > 0, et la trajectoire est une ellipse si 
a = 0. 


3 Systèmes linéaires d’ordre 1 avec second membre 


Donnos-nous une fonction f : R — ÆE continue. Un système différentiel du premier 
degré avec second membre est une équation différentielle de la forme 


(E) z'(t)=ur(t)+f(). 
L’équation sans second membre associée est 
(H) z'(t)=ux(t). 


Elle est encore appelée équation homogène associée. 


3.1 Résultats fondamentaux 


Théorème 8 - Méthode de variation des constantes de Lagrange - 
Il existe au moins une solution de l'équation (E) de la forme x (t) = e“* (c(t)) où c(t) est 
une fonction à déterminer. 


Preuve : On a x’(t) = ue"“t(c(t)) +e“t(c (#)) et 
(E) & ec (0) = F0 & d'( =e *(f (0): 


La fonction # + e-“!(f (t)) est continue, et il suffit de prendre c(t) = fie“ (f(t)) dt 
pour conclure. = 


Théorème 9 Toutes les solutions de (E) sont obtenues en additionnant une solution 
quelconque de l'équation sans second membre (H) à une solution particulière de l’équation 
avec second membre. 


Preuve : Si y(t) est une solution particulière de (E), 


…. a'(t) = ux (t) + f(t) 
a'(t) = ux(t) + f(E) a 


) 
& x — solution de (E}&IXEE x=y+e"t(c).m (x) 


Corollaire 1 -Unicité du problème de Cauchy - 


Il existe une et une seule solution de (E) satisfaisant la condition initiale x (to) = to. 
C’est 


a (4) = y (6) +46 10) (x0 — y (to) . 


Preuve : Les solutions de (ÆE) sont les fonctions x (t) = y + e“{(c) où c € E et où 
y (t) désigne une solution particulière de (E). On détermine le vecteur c de façon à avoir 
x (to) = xo. On trouve ro = y (to) + e"t0 (c) d’où c = €" (xo — y (to)). M 


Le Théorème 9 montre combien la recherche d’une solution particulière de l'équation (E) 
est importante. Dans la pratique, la recherche d’une solution particulière de (E) se fera 
soit par la méthode de variation des constantes de Lagrange (Théorème 8), soit par la 
méthode d'identification exposée dans la Section suivante. 

On remarquera que la recherche d’une solution particulière de x’ (t) = ux (t)+ f (t)+g(t) 
peut s’effectuer en superposant les solutions. Cela signifie que l’on cherche d’abord une 
solution particulière y(t) de x’ (t) = ux(t) + f (t), puis une solution particulière 2 (t) de 
x (t) = ux (t) +g(t), pour finalement constater que y (t) +2 (t) est solution particulière de 
x (t) = ux(t) + f (t) + g(t). Un exemple complet de résolution est donné à la Section 4. 


3.2 Second membre de la forme e“Q (t) 
Théorème 10 Soient u EC et Q(t) € Eft]. L'équation différentielle 
(E) 2'(#) = ux(t) +e#Q(t) 


admet au moins une solution particulière de la forme et P (t) où P est un polynôme de 
degré deg P — deg Q + mult (x), l’entier mult (u) désignant l’ordre de multiplicité de y 
dans le polynôme caractéristique de u. 


Preuve : Notons r = mult (u), Q(t) = > af .et P() = > pet® où p4 et qx désignent 

des vecteurs de E. La fonction x (t) — HP (4 ) est solution d (E) si, et seulement si, 
gel P (+) + et P! (4) = uelt P (t) + e#tQ (+) 

ce qui équivaut successivement à 


(u— uld) P(t) = P'(t) -Q() 


m+r 


J(u— u1d) (pr) t* = D [(& +1) peu — alt 
k=0 k=0 


VRE{O,Lm+r} (und) (px) = (k+1)pr qe (4) 
e Sir —0, l’endomorphisme u — 14 est inversible et (+) équivaut à 


Pm = (u — HId) ” (—Gm) 
Pm-1 — (u 5. ad) (MPm nd dt) 


Dans ce cas, le polynôme P est unique. 


e Sir > 0, décomposons px et gx; dans la somme directe E = N(A\)@...@ NX) des 
sous-espaces caractéristiques : 


{ Pk = Ph + + DR 

Uk = ki + ee + GR, 

Supposons u = À, et notons u; la restriction de u à N'(À;). L’équation (x) équivaut à 
Vi VkE{0,1,.,m+r} (ui —uld)(prx) = (k+1)prrix — dx. (++) 

Si à Z 1, l'endomorphisme u; —114d est inversible et le système (++) possède une solution 


d’après le premier cas. Si à = 1, alors (u1 — u1d)" = 0 et il suffit de choisir po, quelconque 
et définir les py, x, = Pr, tels que 


VR EL LME Hi u1 — Ud) (pu) + ul 


1 
se 
k+1 
pour que les m + 1 premières équations de (+*) soient satisfaites. Ensuite, pm+1,1 étant 
imposé, il suffit de voir que l’on doit seulement vérifier 


Vkefm+l,..;m+r} Prin = (u1 — u1d) (pru) 


1 
k+1 
C’est possible. On remarque surtout que 

1 
puisque (u1 — u1d)" = 0, ce qui prouve que le nombre de coeficients p4, est fini et achève 
la preuve. On a ici trouvé une infinité de polynômes P répondant à la question. mn 


Pm+r+iyu — (ui — md) (Pm+1u) = 0 


Remarque : Le Théorème 10 peut être utilisé dans les cas où le second membre f (t) est 
de la forme Q (t) cost, Q (t)sint, Q(t)cht ou encore Q (t)sht. Il suffit de superposer les 


e “te * 


solutions en notant que cost — etc. 


10 


4 Exemple 


Cherchons les solutions complexes du système différentiel avec second membre 


A — 5x — 3y — 32 
(E) & y = —-3r+5y+3z+e 
MR 9x — Jy — 7z. 
On doit d’abord chercher les solutions générales du système sans second membre 
x — 5x — 3y — 3z 
(H) 4 y = —-3x+5y+3z 
2 = 9x —9y—72. 
Le système (H) s'écrit X’ = MX où 
5 —3 -3 “à x (t) 
M= | -3 5 3 Xl y = Tv 
9 —9 —7 2 z(t) 
Le polynôme caractéristique de M est 
5 — À 3 3 
ON 2 = —X+3X2-4—-(X+1)(Àx — 2). 
9 —9 —7— À 


On sait que toute les solutions de (7) seront de la forme X (t) = Ae-t+(B + Ct)e?t où À, 
B, C sont des vecteurs-colonnes de C?. Pour déterminer les vecteurs À, B, C', on injecte 
cette forme générale des solutions dans l’équation différentielle pour obtenir 

A F4 Ce*.+2(B LCte* M Ae t + (MB + MCt) e?t 

A6 LB LCYE" LICte?t = MAëéTiMBe* + MOCtei 


La famille des fonctions{t + t7e%t} étant libre, on déduit 


meN,acC? 


MA =-A (M+DA =0 (M+I)A  =0 (1) 
MB =2B+C &4 (M-21B =C &4 (M-21B =C (2) 
MO 20 (M=2F)C =0. (M-21Ÿ B =0. (3) 
Il s’agit donc de résoudre (1) et (3), puis de déduire C' par (2). On trouve 
6 —3 —3 zx 0 x —Q] 
(A)s | -3 6 3 wie! 0 Pre A=l 6 = a] où a EC. 
9 —9 —6 2 0 2 —3a] 
De même 
SR 0 _9 9 9 x 0 
aie lé ce, 3 pie lies 9, ag 0 y |=|0 
9 —9 —9 Z 0 —27 27 27 Zz 0 
T bi + bo 
& B—- y = bi où b1,b2 E C 
2 bo 
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puis 


3 —3 —3 bi + bo 0 
C=— —3 3 3 bi L 0 
9 —9 —-9 bo 0 


En conclusion, les solutions du système (H) seront 


—ae t + (b1 + b2) e?t 
X (t) —= Ae t + (B + Ct) et = ae + bie?t 
—3aie-t + be? 


où @a1,01,02 € C. Comme on s’y attendait, l’espace des fonctions solutions de (H) est de 
dimension 3 sur C, et une base de cette espace est 


_et e2t e2t 
re ee s tmlet |; tn 0 
—3e"t 0 ee 


Les solutions générales de (E) s’écrivent comme somme d’une solution particulière de (E) 
et d’une solution générale de (H). Cherchons donc une solution particulière de 


0 
X'= MX + Jet où J = | 1 |. (+) 
0 


Comme —1 est racine d’ordre 1 de Xy (À), on sait que cette solution particulière peut être 
recherchée sous la forme X = (A + Bt)e-t. En remplaçant dans (x), on trouve 


(B=A = ie tr = MÜAEBdE te Jet 
B-A-Bt = MA+J+MBt 
soit 


(M+DA=B-J 
un = 


Ce système équivaut à 


B=(M+IA+J (4) 
Folie. (5) 
6 ES NT be = 0 
5) | -3 6 3 y |=| -3 6 3 21 
9 —9 —6 2 ÿ 0 26 0 
18 —-9 —9 7: —a— À 
& | —-9 18 9 y |=| -6 | & A— a 
EE 2 9 —3a—1 
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Prenons par exemple À — (—i, 0, —1). Alors 


.+# 0 1 
3 
B=(M+I)A+J = 0 +1 = | —-1 
—1 0 3 
Une solution particulière de (E) sera donc donnée par 
44e" 
= (A+ Bt)e = —te t 
(—-1+3t)e t 


et les solutions générales de (Æ) seront 


—aie Ÿ+ (bi +b)eË+(t-5)et 
X= dié Abe rte 
—3aje ? + be? + (3t— 1)e-Ÿ 


On notera que la même méthode peut être utilisée pour résoudre une équation dont le 
second membre est plus compliqué, comme par exemple 


x! — 5x — 3y — 32 + cos 2t 


(E") 4 y = -3r+5y+32+e t 
“Re — 9x — 9Jy — 72. 
Comme cos 2t = HT, le système (E’) s'écrit 
e2t er?t 
; : 2 "2 
XX =MX+Je"+K+L où K = 0 et L — 0 
0 0 


et une solution particulière de (E") sera obtenue en additionnant les solutions particulières 
de chacune des équations X/ = MX +Je !, X'=MX+Ket X'= MX +L. 


5 Equations linéaires d’ordre n à coefficients constants 


Il s’agit d'équations de la forme 
CE) 20) (0) + ana (6) +. + aox (0) = f (D) 


où les a; sont des constantes complexes et où f (t) est une fonctions continue de R dans 
C. Toute fonction n fois dérivable x (t) de R dans C vérifiant (Æ) est appelée solution de 
l’équation différentielle (E). Le second membre de (E) est f(t). L’équation différen- 
tielle homogène associée à (Æ) (encore appelée équation sans second membre) est, 
par définition, l’équation 


(H) a (4) + an 1e 0 (4) + + a0x (t) = 0. 


Définition 1 L'équation À + an_1X" 1 +... + ao = 0 est appelée équation caractéris- 
tique de (H). 
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5.1 Utilisation du système différentiel d’ordre 1 associé 


Théorème 11 1) L'ensemble Sn des solutions de (H) est un espace vectoriel de dimen- 
sion n sur R (ou C). 
2) La fonction x (t) est solution de (H) si, et seulement si, 


P 
æ (8) = SN pi (t)e** 
i=1 


où les À; désignent les racines de multiplicité r; de l’équation caratéristique PURES (ee 
… + ao = 0, et où pit) ERft] (ou C[t]) est un polynôme de degré deg p; < ri. 

3) Le système (tie) icicpocser, C5 une base de SH. 

4) La solution générale de (E) est obtenue en additionnant une solution particulière de 
(E) et la solution générale de l’équation sans second membre. 

5) Il existe une et une seule solutions de (E) satisfaisant les conditions initiales 


(x (to) ,&' (do) ,…, xD) (to) EP 
Preuve : 1) x(t) est solution de (AH) si, et seulement si, 


X (6) = (an (6) , mo (6), En (6) = * (x (6),2° (6), œ0r D) (®) 


est solution du système différentiel d’ordre 1 : 


a (6) 0 1 0 1 (£ 
æ (t) 0 1 2 (t) 
(a - , 
0 1 : 
HA (t) —40 . 7 An-2 —An-1 Tn (t) 
, 
X'() À X(0) 


L'application 
D: SH — Sy 
xt) + X(t)="{xf),zx (©), ne de (#)) 


est un isomorphisme d'espaces vectoriels de de l’espace Sx des solutions de (A) sur l’espace 
Sy des solutions de (H'). Par suite dim Sy = dim Sy = n. 


2) Le polynôme caractéristique x14 de À est 


—À 1 Ji 0 
—À 1 
XA(N =! : : = (1) (NU + an 1 TT +. + 00). 
— À 1 
—&@0 os 7 An-2 —An—-1 — À 
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Cela se montre par récurrence en développant le déterminant suivant la première colonne : 


—X 1 0 
—} 1 
XAU) = (à)! : , — ao(—-1)"" 
— À 1 
—@1 F0 —An—2 —An—1 — À 


= (=) (AT Han NT +. + «1) + ao (—1)" (hypothèse récurrente) 
(—-1)" (x + HAT ne ao) ; 


p 
On sait que la forme générale des solutions de (H') est X (t) = 3° P, (t) e\it où 
i=1 
P,(t) € Eft] et où deg P; < r;. Compte tenu de l'identification du 1), x (t) aura bien la 
p 
forme x (t) = S pi (t) et où pi (t) ER[t] (ou C[t]) et deg ps < ri. 
i=l 


3) Le système de n vecteurs (ésiehit) engendre l’espace vectoriel Sy des solu- 


1<i<p,0<s;<ri 
tions qui est de dimension n d’après 1). C’est donc une base de Sx. 

4) Même démonstration qu’au Théorème 9. 

5) Fixer les conditions initiales (x (40) ,& (to), 25200 (to)) = (æ1, T2, .…., Tn-1) revient 
à imposer la condition initiale X (to) = (x1,%2,...,%n-1) à la solution de l’équation du 
premier ordre X’(t) = AX (t)+F(t), où F(t) = t(0,0,..., f (t)). Le résultat provient de 
l’unicité du problème de Cauchy pour l’équation X’(t) = AX (t) + F(t) (Théorème 1). 


5.2 Autre méthode 


Si l’on sait que l’ensemble Sx des solutions de l’équation différentielle 
(H) at (4) + an_12 0) (4) +. + aox (€) = 0 


est un espace vectoriel de dimension n, on peut montrer directement le résultat suivant 
(en conservant les notations de la section précédente) : 


Théorème 12 Une base de l’espace Sx est B — (ésieit) need 


Preuve : dim Sy = n = #B de sorte qu’il suffise de montrer que B est un système libre 
et que chaque t‘ie\it est solution de (H) 

e Montrons que t‘e\? est bien une solution de (H). Pour cela, supposons que À soit 
racine de multiplicité r de l’équation caractéristique, et que s < r. Soit L le polynôme 
L(1) = + ap_1 À 1 +. + ap. On pose a, = 1 et l’on définit l'opérateur différentiel 


L’équation (H) s’écrit alors L (2) x (t) = 0 et tout revient à montrer le Lemme suivant : 
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Lemme 3 Sip(t) est un polynôme et si À € K, alors 


1(5) GO) = 52020 (0 


Preuve du Lemme : 
L (à) (» (t) ex) = Du (» (t) 4 AS 43 cn D) ( px 


(k) il (&) n k 
x P D sh à p® (4) de Là 
_ DD LE EL = k! D_ aie (6) la 


k=0 i=k k=0 i=k 
mer (6) Cr LL MNT PC) ,&) 
= ASE ( GE «) h=e re DC 
k=0 i= k=0 
e Le fait que B est un système libre provient du 
Lemme 4 Si À1,.…., À, sont n complexes distincts deux à deux et si p1 (t) ,.…., Pn (t) sont 
des polynômes de K [t], alors 
Ait Ant. 2 _ _ _ 
pmitt)e +... +m (te Spies = (00: 


Preuve du Lemme : Raïsonnons par récurrence sur n. C’est trivial si n = 1. Au rang 
n, l'équation p1 (t) et +... + p, (t) ent = 0 devient 


p1 (t) cb pet (t) eln-1t + p, ()=0 


où les coefficients y; = À; — À, sont non nuls et distincts deux à deux. Si p, (t) est nul, 
il suffit d'appliquer l’hypothèse récurrente pour conclure. Sinon posons d — deg p, (t) et 
dérivons l'égalité d + 1 fois. On obtient 


n—1 
> (pi () etti 70 0 
i=1 
soit 
n—1 d+1 e 
D 2 Chan (£) DEEE — 0 
i=1 k=0 


L'hypothèse récurrente entraîne pC pi ) GE) pÉTITE 2 G pour tout 4. Si p; (#) n'était 
k= 


pas le polynôme nul, alors le ee (pi (6) , pt (t), … plie? 9 (t)) serait libre puisque 
formé de polynômes de degrés étagés, et l’égalité précédente entraînerait HET = (0, soit 
1; = 0, ce qui est absurde. On à donc p; (t) = 0 pour tout à € N,,_1, et cela entraîne aussi 


Pn (t) = 0. La récurrence aboutit. m 
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5.3 Second membre de la forme el“#q(t) 
Théorème 13 Soient q(t) un polynôme de Cft] et u € C. L'’équation différentielle 
(E) 20 (4) +an12 070 (6) ++ aox (t) = q(t) e4 


admet une unique solution de la forme x + p(t) e“t où p(t) est un polynôme de valuation 
val (p) = mult (4) et de degré deg p = deg qg + mult (4), où mult (4) désigne la multiplicité 
de u dans le polynôme caractéristique de (E). 


Preuve : Posons d = deg q, s = mult (y) et q(t) = qo + … + qat®. D’après le Lemme 3, il 


d+s 
s’agit de trouver un polynôme p(t) = pstf +. + parstTts = ST pytf tel que 


k=s 

d LR MONET 

— nt) — Lht 2222) (= ut 
L (7) Ge) DE =abet 

Comme LÜ) (x) = 0 pour tout à < s, cela équivaut à 

n dt ; 
ss "0 (E—1)..(k—i+ DE = (6) 
i=s k—i 


D GED Gr Pdf) 


i=s j—=0 
d |inf(d+s,n) - 
LG) De. 
de Es G+pG+5-0.G+ 016 200. 
3=0 i=s 


En fait, dans cette dernière somme, 4 varie de s à d+s. En effet, p étant de degré < d+5, 
on a p;+; = 0 dès que i+3j > d+s, donc à fortiori si à > d+s. Il faut maintenant résoudre 


d+s 
: 5 nn : 
Vj € {0,.…,d} pis (+ D + 1). (+ D =. 
il 


i=s 


Ces d équations forment un système triangulaire dont les inconnues sont les ps .…., pa+s et 
dont les coefficients des termes de la diagonale principale sont non nuls (en effet L(®) (4) Z 0 
puisque x est de multiplicité 4 dans le polynôme L(t)). Ce système est inversible est 
détermine les ps.….,pa+, de façon unique. 
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